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The unitary transformation that leads from the minimal-couphng description to the electric- 
dipole one is analysed in detail. The momentum cut-off function f{k), which is understood in the 
definition of such a transformation, is obtained explicitly by use of the variational method. We 
give an expression for f{k) in terms of the electron-density and of the mean-squared value of the 
electric-dipole moment of the atom in its ground state. A relevant consequence concerns the van 
der Waals interaction, whose long-distance behaviour turns out to be proportional to R~^. 
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I. INTRODUCTION 

In non-relativistic quantum mechanics the standard 
description for an atom interacting with the radiation 
field is provided by the minimal-coupling (MC) Hamilto- 
nian. As it is well known, the MC Hamiltonian can be 
obtained in a natural way as a non-relativistic limit for 
the Dirac Hamiltonian (see e.g. Ref. [l| ). Nevertheless 
it is known also that, in the long- wavelength limit, the 
gauge coupling of the MC description is replaced usually 
by the electric-dipole (ED) coupling — d • E. 

The transition from one description to the other one is 
analysed, in several books, in terms of classical electro- 
magnetic field (see e.g. Ref. [1] [1] ). Much more elegant 
is a method in which, from the beginning, the electro- 
magnetic field is introduced as a quantum field (see e.g. 
Ref. 4]). Such a procedure is based on the introduction 
of a unitary transformation, which leads from the MC 
representation to the ED one. Some interesting features 
concerning such a transformation will be analysed in this 
paper. 

Let us introduce the problem briefly. The total Hamil- 
tonian Ht, which describes an atom interacting with the 
electromagnetic field, consists of 



Ht — Hp 



H'm- 



(1.1) 



In eq.(l.l) the free- field Hamiltonian Hp is given, in 
Gaussian units and in the Coulomb gauge, by 



Hf j dr:{El+B^) 



(1.2) 



(semicolons mean normal ordering). The Hamiltonian 
H'j^, which describes both the atom and its MC interac- 
tion with the radiation field, has the well-known form 



H',, = ^■.\-^hV + -A{r)\^ 
Zni c 



-V{r-R). (1.3) 



For the sake of simplicity we assume a single electron, 
we neglect the spin and we assume that the centre-of- 
mass of the atom be fixed at the position R. 

In eq.(1.3) the normal ordering, which concerns the 
term only, has the effect of fixing a reference-level for the 
energy, according to the following relation 

(5|(0|i/;,|0)|5> = {g\HMo\9) = Eo. (1-4) 
In eq.(1.4) |0) is the vacuum state for the field, 
ft2 



HMo^-^y^ + V{r-R) 



(1.5) 



is the Hamiltonian for the unperturbed atom and Eq is 
the energy-eigenvalue corresponding to the ground-state 

The transition to the ED representation is accom- 
plished through the following unitary transformation [3| 



U = exp[- — d- A(H)] 

he 



(1.6) 
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where a high-momentum cut-off is understood implicitly. 

In this paper we will proceed in the following way. First 
of all we will introduce a cut-off function f{k), for the 
photon momentum k, explicitly in the definition of the 
unitary transformation t/ ofeq.(1.6). This will be done in 
eq.s (2.1) and (2.3). Subsequently the function /(fc)will 
be determined by the variational method. As a result we 
will obtain for f{k) an expression involving the electron- 
density ne{k) and the mean-squared value {(P) of the 
ED moment of the atom (see eq.(3.8) in the sequel). In 
the ED representation the interaction-Hamiltonian will 
contain the cut-off function /(fc) explicitly and it will 
differ significantly from the simple expression —d-E±{R), 
which is used currently. 

The case of N atoms is treated in Sec. IV . The uni- 
tary transformation is redefined suitably in eq.(4.1). In 
this case the transformation gives rise to a dipole-dipole 
interaction among different atoms, to be ascribed to the 
transverse field. Such an interaction overlaps with the 
dipole-dipole interaction due to the Coulomb field. As 
a result, at large distance an exact cancellation of the 
contributions proportional to occurs. So the leading 



contribution to the large-distance dipole-dipole interac- 
tion turns out to be proportional to R^^. This would 
imply for the van der Waals potential a leading contribu- 
tion proportional to and not to R~^. This is shown 
in Sec.V. Actually the experimental status does not look 
conclusive. 



II. THE ELECTRIC-DIPOLE 
REPRESENTATION 

The transition to the ED representation is achieved 
through the following unitary transformation (Gaussian 
units are used) 



where 



U = exp[-—d- Af(R)], 

he 



d = -e{r - R) 



(2.1) 



(2.2) 



is the ED operator. In eq.(2.1) a momentum cut-off has 
been introduced, by defining 



Af{r) = cJ^Vffl[e,<,a.„e^'='- (2.3) 
In a similar way we define 



^{(0 = *\/^E/>)v^tefcaafe„e*'=-'- (2.4) 



V 



ka 



^koiO-ka^ J- 

The true field operators A and E±_ correspond to / = 1. 

The real isotropic function /(fc) carries out the high- 
momentum cut-off. It will be determined in Sec. Ill by 
the variational method. 

The presence of the cut-off function f{k) affects the 
commutation relations among field operators, namely 



[Al{R),E^,{r)] 



(2.5) 



kikj , 



Owing to this result, the transverse electric field oper- 
ator transforms as 



An 



(2.6) 



ki 



The vector-potential A, as well as the magnetic- 
induction B, is not modified by the transformation. 



These results entail that, in the ED representation, the 
free-field Hamiltonian Hp reads 

UHpU-^ = Hp (2.7) 

-Y.d^y-'Y. J drE^j{r)e^'^<^-^^m{Sij - ^) 

i ki 

+27r^d.d,y-^^f(fc)(5.,-^). 



By calculating the integral in the second term we obtain 
E / drE^j{r)e^''<^-^)fik){S,j - ^) (2.8) 



= I 



2TTh 



J2 hk)y/^[ekccakae^^-^ - el^a+^e 



ha 



(2.9) 



= eUr)- 

The integral in the third term gives 
2.y-iE/'W(^^i-^) 

k 



= ^^-jj^ dkk'fik). 

In performing the calculation, the scalar contribution 
{L = 0) has been separated from the null-trace one 

{L = 2). 

In this way we obtain 

UHpU-^ (2.10) 

= Hp -d-EUR) + — / dkk^P{k). 
371" Jo 

Finally, putting all these results together, we obtain 



UHtU-^ 



(2.11) 



2m 



ih"^ + ^[A{r) - A^ {R) 



+V{r-R)-d-E{{R) + —\d\^ dkk^f{k) 



+^(oii^^(H)no) 



:{0\A{r)-AfiR)\0). 



We have to explicitate the last two terms. From the 
first one we obtain 



(oiiA/(i2)no) 



(2.12) 



= 2-Ke^KV-^ y k-'fik) ~ ^ / dk kf{k), 
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where Ac = h/m,c is the Compton wave-length for the 
electron. From the second term we obtain 



;{Q\A{r)-Af{Rm (2.13) 



{r-R) 



(2.14) 



In this way from eq.(2.11) we obtain 

= HF + ^:\-mV + -[A{r) - Af{R)]\^ : 
2m c 

+V{r - R) - d ■ e{{R} + —\df / dkk^P{k) 



ik-(r-R) 



+^ / dkkf{k)-^ / dkk-'f{k)e 



Currently eq.(2.14) is written with / = 1, and the last 
two terms (the first of which would diverge for / = 1, 
while the second one would give the result — 2e^Ac7r~^ 
\r — -R|~^) are neglected. In the \d\^ term, which di- 
verges also for / = 1, a momentum cut-off is understood. 
Nevertheless eq.(2.14) differs from the expression 



UHtU-^ =Hf + Hmo - d ■ E^{R), 



(2.15) 



which is assumed usually. 

In the sequel we will refer to eq.(2.14), where the func- 
tion / will be determined in Sec.III by the variational 

method. 

Let us proceed to a preliminary analysis of eq.(2.14). 
By averaging both sides of eq.(2.14) over the vacuum 
state of the field we obtain 



H'mo = {0\UHtU-'\0) 



(2.16) 



9 /"^ 

e2 



2 A, fOC 

/ dkkPik) 

dk k^Pik) 



27r2 



j dkk-^f{k)t 



ik-{r-R) 



This equation can be generalised easily to the case of a 
many-electron atom. 

The hydrogen- atom case is simple. In fact in such a 
case eq.(2.16) reads 

H'MQ = HMQ + V{\r-R\), (2.17) 

where the additional potential V has the form 

V = Co + \mnl\r - + V"(|r - R\). (2.18) 



In eq.(2.18) Co is the following constant 



Co = TT-h'^Xc dkkf{k), (2.19) 
Jo 

while, for the harmonic-oscillator potential in the second 
term, one has 

1 o 2 poo 

-mQl = — dkk'^fik). (2.20) 

2 3n Jo 

Finally the last term consists of the following additional 
potential 



V"(r) 



277 

2e2A, 



nr 



j dkk-^f{k) 



dk f{k) sin kr 



e--' (2.21) 







where the last side represents the asymptotic expression 
for large r. In the calculation of the last expression the re- 
sult /(O) = 1 has been used. This result will be obtained 
in the sequel (see eq.(3.9)). We observe that V"(r) does 
not vanish for r = 0. In fact it is 

V"(0) = -27r-^e^Ac / dkkf{k). (2.22) 
Jo 



In the case of a Z-electron atom, eq.(2.2) is replaced 



by 



d=-eJ2ifA-R) 



(2.23) 



and eq.(2.14) reads 

UHtU-^ (2.24) 

= + i E ^ I - ^^^^ + ;[^(^^) - ■ 



A<B 

-d-E{{R) + —\df dkk^P{k) 
Stt Jo 

POO 

+7r-^e2AcZ / dkkp{k) 
Jo 

I dkk-'f{k)e'''<'-^-^\ 



e^Ae 
" 27r2 



III. VARIATIONAL CALCULATION 

We will determine the function / by the variational 
method, assuming a trial state of the form 



|Vt) = C/|0)|<?), 



(3.1) 
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where \g) is the ground-state for the unperturbed atom. 
The corresponding wave-function will be indicated as 



Mr-R) = {r\g). 



(3.2) 



We observe that the trial state of eq.(3.1) is normal- 
ized. For the sake of simplicity we assume that the wave- 
function <j)o is isotropic {I = 0) and we neglect the spin. 
The energy of the system, averaged over the state (3.1), 

£o[,/] = (V'Tli^Tl^T) = {gmU-^HTU\Q)\g), (3.3) 

depends on / functionally. Using eq.(2.16) we can write 

(3.4) 



^o[/] = So +7r-^e^{Ae / dkkf\k) 
Jo 



r dkep{k)-2K 

"J e Jo 
In eq.(3.3) {(]?) is given by 



dk kne{k)f{k)}. 



{d^) = e^ J drr^\Mr)f 

poo 

= 47re2/ drr^|</)o(r)|2, 
Jo 



(3.5) 



while 



ne(fc) = J dre'''-^\Mr)f 

/>CXD 

= 47rA;~^ / dr r|(/)o('')P sinfcr ~ 1 
Jo 



(3.6) 



6e2 



consists of the Fourier-transform of the electron-density 
ne{r) = \4>o{r)\'^. The last side represents the asymptotic 
behavior for A: — >■ 0. 

Eq.(3.4) confirms the need for a high-momentum cut- 
off. In fact it can be seen immediately that for / = 1 the 
energy £q diverges. 

The function f{k) can be obtained by minimising the 
energy Sq of eq.(3.4). We observe that a vanishing / 
would return the result Sq = Eq, in agreement with 
eq.(1.4). Therefore it can be expected that any non- 
trivial solution / to the variational problem would corre- 
spond to an energy £q < Eq. 

By minimising the energy £o[f] of eq.(3.4) with respect 
to f{k) we obtain the equation 



= d£o[f]/Sf{k) 

= -ke^iiXc + ^k)f{k) - XcMk)]. 
The solution to this equation is given by 



(3.7) 



/(fc)=ne(fc)(l + 



3e2A< 



k)-\ 



(3.8) 



It can be seen immediately that /(O) = 1, as it could 
be expected. Furthermore the following asymptotic be- 
haviour for small k 



fik) 



(3.9) 



can be derived immediately from cq.s (3.6) and (3.8). 
Moreover from eq.s (3.8) and (3.6) it follows that, for 
fc ^ oo, f{k) vanishes as k~^. 

Finally, using the result of eq.(3.8) in eq.(3.4), we ob- 
tain for the energy £q the result 

/■OO 

£o[f] = Eo- TT-^e^Ac / dkkf{k)he{k) (3.10) 
Jo 

= Eo- TT-h^Xe dk khl{k){l + |^fc)-\ 

with / given by cq.(3.8). Eq.(3.10) shows that £o[f] < 
Eq, as expected. 

We observe that the second (negative) contribution in 
eq.(3.10) would have the effect of shifting every energy- 
level of the atom by the same amount. Therefore such 
an effect cannot be observed through any spectral-line 
analysis. 

As an example, for the hydrogen atom in the Is state 
we have 



f>o{r) = (7ra^)-^e 



3\-lp-r/a 



an correspondingly 



e-'{d') = 3a^ 



(3.11) 



(3.12) 



where a = H'^/me'^ is the Bohr radius. Putting = ^ka, 
we have 

he{k) = (1 + zl)-^ ~ 1 - ^fc^a^ + ... (3.13) 
In this case from eq.(3.8) we obtain 

f{k) = {l + zlr'{l+Aa-^z,r^ (3.14) 
= 4a-\l + zl)-^{zk+a/4)-\ 

where a = /he is the fine-structure constant. 

Let us evaluate how many photons are present in the 
state (3.1). Prom eq.(2.1) one obtains 



UttkaU ^ = aka + i 



2n 



e-"=-«/(fc)d.eL (3.15) 



hckV 

and a similar expression for a^^. In this way we obtain 

= J2(9\maL^kaU-'mg) (3.16) 
27r 



hckV 



f'{k)[{d')-k-'{{d-kr)]. 



In the continuus limit, by integration over the direction 
of the wave- vector k, one obtains 

nph{k) = ^ fdih k^nf - ^kp{k), (3.17) 
Stt J 6-Knc 

with a total number of photons given by 



ph 



POO 

j dkuphik) 
Jo 



'Snhc 



dkkfik). (3.18) 
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For a hydrogen atom eq.(4.12) gives the result 
1.4 X 10"^. 



(3.19) 



For the calculation of the integral see the Appendix A. 

It is interesting to note that the maximum value for 
Uph is attained for z = zq = 1/7, e.g. for A ~ 22a, in the 
far-ultraviolet region. 

The transformation U of cq.(2.1) could be introduced 
in the frcc-clcctron case also, in order to seek for a pos- 
sible spontaneous localisation for the wave-packet, due 
to the second term in eq.(3.10). For e.g. a Gaussian 

(^0 e-'^/''\ cq.(3.10) would give Ai;o l^i^f- 

But for a non-relativistic electron this negative contribu- 
tion to the total energy would be too small to have any 
observable effect. 



IV. N ATOMS 

In the case of N identical atoms, the unitary transfor- 
mation U of eq.(2.1) becomes 



C/ = exp[-^^da-A/(i?„)]. 



N 



(4.1) 



a=l 



As a consequence, eq.(2.6) now reads 



[UE^iU-\r) = EUr 
An 



(4.2) 



o 3 



The Hamiltonian Hm of eq.(1.3) is replaced by 

^M=V[;^:|-i?iVa + -A(ra)|2: (4.3) 

a 

+V{r - Ra)] + 5^ ^ rfai QijiR'"'') dbj. 

a<b ij 

In eq.(4.3) we have put = Ra — Rb. Furthermore 
the matrix g is given by the expression 



g,,iR)=R-^iS,j-3^: 



(4.4) 



An 
'V 



A;2 



The last expression is correct for R> 0. 

In eq.(4.3) the last term represents an interaction, due 
to the Coulomb field (longitudinal electric field), among 
the ED moments belonging to different atoms. In the 
standard treatment, this term alone is responsible for 
the van der Waals interaction. In second-order perturba- 
tion theory it gives rise to an interaction-potential, whose 
leading contribution at large distance is proportional to 



R~^. However we will show that the unitary transforma- 
tion of eq.(5.1) gives rise to a similar interaction term, to 
be ascribed to the transverse field. In fact we have 

UHtU-^ = Hf (4.5) 

a 

+V{r, - Ra) - d„ • EiiR,) + —\da\^ dkk^m 



+7r- 



Jo 



dkkp{k)} 



+ EE^-ri,(ii"'')4,. 

a<b ij 

The last term represents the interaction among the ED 
moments referring to different atoms, due both to the 
Coulomb field and to the modified transverse electric- field 
E^^. The matrix F is given by 



(4.6) 



, A-K 



{5,,V^ - d,di)-Y,k-^[l- f\k)V 



k R 



2 /""^ 

^ {dijV^ - didj)- dk[l-f{k)] 

7^ Jo 



sin kR 
kR 

The leading contribution to this expression can be eval- 
uated by use, in the calculation of the last integral, of the 
asymptotic expansion for /, given in eq.(3.9). We obtain 



8{d^) „ 2 P , . ^d^) r. 2 



(4.7) 



where the two contributions proportional to R~^ cancel 
each other. In this way, for the leading contribution to 
Tij we obtain the expression 

rM-i^R-H^^-s.,). (4.8) 



37re^Ac 



i?2 



For example in the case of a pair of hydrogen atoms 
eq.(4.8) reads 



T,j{R) 



'%(|)^(4^-^.)- (4.9) 



iraa 



The interaction due to the last term in eq.(4.5) is re- 
sponsible for the van der Waals interaction. This will be 
analysed in the next section. 
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Quite relevant in the A''-atom case is the fact that the 
photon number of eq.(3.18) is multiphed by N. In fact a 
uniform distribution of N atoms in a volume V can be as- 
sumed as a rough representation for a refractive medium. 
In such a case, the number of photons per unit- volume 
would be given by 



V 



VSnhc 



dkkfik), 



(4.10) 



with the following distribution with respect to the fre- 
quency UJ 



nphiiu) N {(f) 



V 



V 6nnc-^ c 



(4.11) 



Assuming N/V 10' 



23 



cm 



and using the result of 
eq.(3.19) one obtains Nph/V 10^'^ cm"^. This resuh 
(very roughly) represents the total number of photons 
per cubed centimeter, which are present actually in a 
refractive medium in its ground state. 



VAN DER WAALS POTENTIAL 



Let us recall that both U and U' are negative. From 
eq.(5.4) one can see that \U{R)\ > \U'iR)\JoT R < 9Ax 
lO^a. Actually the experimental status does not look 
conclusive. 



VI. CONCLUSION 

We have shown that, in the ED representation, the 
correct Hamiltonian describing an atom interacting with 
the radiation field is given by eq.(2.24), with / given 
by eq.(3.8). Such an expression appears sensibly more 
complicated than the simple expression —d ■ E±, which 
is assumed usually. 

We have analysed some consequences concerning 
namely the photon-distribution and the van der Waals 
interaction. In the last case we have obtained the re- 
markable result that, at large distance, the van der Waals 
potential between two atoms is proportional to R~^. We 
are not aware of any direct measurement of the van der 
Waals potential, between two neutral polarizable atoms, 
expecially at large distance. A wide bibliography on the 
van der Waals interaction can be found in Ref. [6|. 



The long-distance interaction between two neutral 
atom is due essentially to the van der Waals force. Up to 
the second order the perturbation theory gives an explicit 
expression for the dipole-dipole interaction energy.The 
standard calculation is based on the Coulomb interac- 
tion as given by the last term of ec^(4.3), with gij given 
in cq.(4.4). The well-known result Q is 



U\R) 



(5.1) 



6lf)2 



2Eq — — Eb^ 



which is proportional to R~^. 

However in Sec. IV we have shown that the complete 
dipole-dipole interaction, is represented by the last term 
in eq.(4.5), where the matrix gij is replaced by of 
eq.(4.8) (or by eq.(4.9) for a pair of hydrogen atoms). As 
a consequence, the potential U' of eq.(5.1) is replaced by 



16(^2 



E 

6162 



i?8 37re2Ac 



(5.2) 



2Eo — Ebi — E\ 



f)2 



which is proportional to i? ^ . 

The ratio of the two expressions is 



6i?2 ^37re2Ac' 
For two hydrogen atoms eq.(5.3) reads 



u{R) 

U'{R) 



6 ^R' ^ira' 



B.9 X 10^(-)2. 
R 



(5.3) 



(5.4) 



Appendix A: Integrals 

We outline the method used to calculate the integrals 
referring to the hydrogen atom. 

The integrals are either of the following type 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 



In{A)^ / dz{z^ + Ar\ 
Jo 

or of the following one 

/■OO 

Jo 

with n, m positive integers. 
For n > I we have 

I,,iA)^{~r~'[{n~iy.]-'^^h{A), 



with 



hiA) 



dA" 



dz{z^ + Ay^ = 

2VA 



In a similar way, for n > 1 and m > 1 we have 

Inrn{A,B) 



(n- l)!(m- 1)! 9^"-! dB^ 



-hi{A,B) 



The integral In can be calculated by standard tech- 
niques. The result is 

/•OO 

In{A,B)= dz {z^ + A)-\z + B)-^ (A6) 
Jo 



{A + B 



2n-1 



ttB 
'2v^ 



In A -In B]. 
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At the end of the calculations one has to put A = 1 and B = a/ 4. 
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